Radiation from accelerated impurities in Bose-Einstein condensate 
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We investigate radiation spectra arising from accelerated point-like impurities in the homogeneous 
Bose-Einstein condensate. A general formula for the radiation spectrum is obtained in the integral 
form as a function of given impurity trajectory. The Planckian spectrum is obtained for a special 
accelerated motion, which is shown to be unphysical. Non-Planckian spectrum is found in the case 
of a uniformly accelerated impurity. We compare our result with similar settings as discussed in 
other quantum many-body systems. 
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I. INTRODUCTION 

It is now a well-established experimental fact that con- 
densation of mesoscopic number of atoms can be created 
and manipulated for relatively long time-scale on the or- 
der of seconds. These experimental progress enable us to 
probe properties of ultra-cold gases and Bose-Einstein 
condensate (BEC) from various different perspectives. 
Among them, many authors have been investigating the 
response properties of BECs when time-dependent per- 
turbations are applied [1-7]. On general grounds, time- 
dependent disturbance on quantum many-body system 
will give rise to energy dissipation and hence to destroy 
the order and coherence in the system after some time 
later. However, due to the phenomenon of superfluid- 
ity of BECs, energy dissipation may or may not take 
place when time-dependent perturbation is applied to 
the condensates. An important question is to know how 
and when such a dissipation process occurs. This paper 
studies the radiation spectra accompanied by accelerated 
motions of an impurity in the homogeneous condensate. 
The motions of classical impurity is a special class of 
time-dependent disturbance and can be realized by using 
a detuned laser, a small number of different species of 
atoms, or many other microscopic objects in BECs. We 
are interested in the form of spectra obtained for various 
kinds of accelerated motions of the impurity. 

Recently, the idea of simulating relativistic and cosmo- 
logical effects in non-relativistic quantum (many-body) 
systems has gained great attention from various fields of 
researchers. There have been many activities devoted to 
explore this idea and how to realize it in laboratories. 
Several authors have studied radiation spectra arising 
from accelerated disturbance and responses of acceler- 
ated detectors in several different quantum many-body 
systems [9, 13-16]. They found the Planckian distribu- 
tion for the radiation spectra and for the detector's exci- 
tation probabilities. These results were then interpreted 
with a simulation of the Fulling-Davies-Unruh (FDU) 
effect in which any accelerated observer perceives the 
empty Minkowski vacuum as a finite temperature equi- 
librium state [17-19]. To be more precise, these studies 
concern simulating the DeWitt-Unruh (DU) effect and 
the FDU radiation. The former states that a particular 



accelerated detector is excited in the same way as being 
replaced in a thermal bath at rest and the latter concerns 
the phenomenon of thermal radiation from an accelerated 
particle in the vacuum [20, 21]. Regarding the DU effect, 
there exist several non-thermal response for the detector 
[22-24]. There also exists a debate if the original deriva- 
tion of the FDU effect is mathematically correct or not 
[25] . It is questioned in this paper whether or not such an 
appearance of Planckian distribution in an analog situa- 
tion can be related to the FDU effect. We briefly discuss 
the consequence of our result for proposed emergent the- 
ories of gravity in quantum many-body systems. 



In this paper we investigate a coupling of a classi- 
cal point-like impurity to the homogeneous condensate 
at zero temperature and study the radiation spectra ac- 
companied by the accelerated motions of impurity. The 
impurity is assumed to move along a given trajectory 
without any back reactions. When the total number of 
bosons in the system is very large, the effect of point- 
like impurity can be represented by the time-dependent 
driven oscillator model, which can be diagonalized by 
a suitable unitary transformation. The radiation spec- 
trum is obtained in the integral form as a function of a 
given impurity trajectory. In this paper, the Planckian 
spectrum is also obtained when the impurity is acceler- 
ated in a special manner inside the homogeneous conden- 
sate. However, this special class of impurity trajectory 
is shown to be unphysical and hence would not be real- 
ized to simulate such a black-body radiation. A radiation 
spectrum is also obtained analytically for a uniform ac- 
celeration and they are analyzed in detail. 



The outline of this paper is as follows. Sec. II provides 
a model and its diagonalization within the Bogoliubov 
approximation. We obtain the radiation spectra for given 
impurity trajectories in Sec. III. Two accelerated motions 
are studied and are analyzed with their corresponding 
radiation spectra in Sec. IV. We discuss our result and 
its connection to the FDU effect in Sec. V. 
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II. MODEL HAMILTONIAN 

The model Hamiltonian for weakly interacting bosons 
of mass M is 



H = Jd 3 x %fl(x) 



h 2 v 2 

2M 



+ ■ 



■ J d s x^{x)^(x)i){x)i>{x), (1) 



where the general two-body interaction term between 
bosons is approximated by a contact interaction g8(x — 
x') in the low energy region. The field operators tft (x) — 
ijr(x, t) and ip(x) = ip(x, t) in the Heisenberg picture sat- 
isfy the equal time canonical commutation relations. We 
assume that the repulsive interaction g > and the con- 
dition of diluteness na 3 <C 1, where n is the number den- 
sity of massive bosons. We consider the following local 
interaction Hamiltonian between bosons and a point-like 
classical impurity moving along a given trajectory 



Hi(t) = A / d 3 x S(x - C(t))^(x)ip(x), 



(2) 



where A represents the coupling constant between bosons 
and the impurity. As we discuss later, A can be time- 
dependent in general when the impurity coupling is 
switched on and off in experiments. For the homoge- 
nous bose gas under consideration, the field operators 
are expanded in terms of the plane wave basis with peri- 
odic boundary conditions in a finite size box V = L 3 to 
obtain the total Hamiltonian H = Hq + Hi as 



h = eka i ak + or? 2 a i+ q a i'- q a * a k' 



k,k' ,q 

A 



+ ££ e - ife - cw 4 +fe ^ (3) 



where e k = h 2 k 2 /2M is the free kinetic energy of bosons 
and the summation is taken over integers rij for k = 
2ir(n x ,n y ,n z )/L. In the following we are interested in 
the effects of impurities to the homogeneous condensate 
in the large N limit. We will investigate the first order 
correction on the homogeneous condensate. 

We follow Bogoliubov's treatment to simplify the full 
Hamiltonian within the number- fixed framework [28]. 
Define the operator j3o = (No+l) -1 / 2 ^ where No = a^ao 
is the number operator of condensed bosons. This op- 
erator and its hermite conjugate satisfy (3oPo — 1 an d 
4/3 = 1 - n , where n = \N = 0)(7V = 0| is the 
projection operator onto the state N = 0. In the pres- 
ence of the homogeneous BEC for the zero mode, we 
can exclude the state No — to approximate /?J/?o — 1 ; 
i.e., /3o is an almost unitary operator, and [A),/^] — 0. 
More precisely this approximation holds between matrix 
elements when we compute expectation values. Then 



bare zero mode operators a and ao- We then introduce 
a new set of the creation and the annihilation operators 
aj[, and ctk for non-zero momentum: a\. — a^Po, ak — 
0oO,k for k 7^ 0. They satisfy the equal time canon- 
ical commutation relations: [ctk(t), oi\,i{t)] — Sk,k' — 
avajjlo ~ S k ,k', [ak(t),a k '(t)} = [a' k (t), a\, (t)] = 0. 
The vacuum state of the Fock space Ho for the new set of 
the creation and the annihilation operators is the same 
as the plane wave basis, i.e., afe|0) = 0. 

The total number operator N = N + ^ a\ak com- 
mutes with the total Hamiltonian, where the prime over 
the summation is used to indicate the omission of the zero 
mode for the summation. Therefore we can choose the 
basis in which N is diagonalized with given value N. This 

gives the condition No + u\.Uk = N which allows us 
to eliminate the zero mode occupation number operator 
No- Using an expansion for the square root function we 
can proceed to do a systematic expansion in 1/ \f~N . The 
Bogoliubov transformation to diagonalize the Hamilto- 
nian without an impurity term is Ub = exp(iGs); 



G B = \Y1 6 k a W-k + h - c -> 



% =tanh- 1 (- 



gn 



1 fojJk + £fc + gn 
hw k = y/e k (e k + 2gn). 



(4) 
(5) 
(6) 



Bogoliubov's excitation (the "bogolon") is created and 
annihilated by the operators 

b k = u Ba\U B = a^cosh^ fe + a_ fc sinh6» fe , (7) 
bk = U B akU B = a k cosh9 k + cJ_ k sinh6> fe , (8) 

respectively. Keeping up to the second order in the new 
set of creation and annihilation operators, we obtain the 
effective Hamiltonian for the impurity motion in the ho- 
mogeneous BEC as 

H = Eo+Y! f^kb{b k + (fk(t)b{ + h.c). (9) 

Here E = gnN/2 — gn/2 + An + \ (haj k — e k — gn)/2 
is the ground state energy without impurity and 



fk(t) = nX 



Nhcok 



ife-C(t) 



(10) 



fNo0Q and a — PoVNq allow us to eliminate the 



The homogeneous condensate is defined by the new vac- 
uum state of the Fock space Hb- b k \BEC) = 0. 

Thus, the impurity effects in our model are described 
by a time-dependent driven oscillator model, or the time- 
dependent van Hove model. It is noted that the effec- 
tive interaction Hamiltonian, the third term in eq. (9), 
is changed from the density coupling (2) to the linear 
coupling; \fn\ J dx 3 S(x — <^(t))(tp^ + -0), where tp is the 
bosonic filed operator excluding zero mode. Note that 
the factor 1/y/N in Eq. (10) is due to the condensation 
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which introduces the significant difference and gives rise 
to the macroscopic effect. 

We can solve the above Hamiltonian (9) exactly with 
the initial condition such that the operators evolve with- 
out impurity asymptotically at the infinitely remote past 
t = — oo. The solution is 



b k {t) = V£{t) + 4> k {t) 

where <j>k{t) is a c- number function 



m) + 



inX 



h Y NhjJk 



Ik(t)e 



-iuikt 



h(t) = f 

J —I 



jj£ e iuj k t' -ik-C(t') 



(ii) 

(12) 



(13) 
(14) 



When the impurity is injected at some time ti, the solu- 
tion is given instead by the similar formulae (11,12,13,14) 
with the replacement of the time integral /_ — > J t .- 

The asymptotic creation and annihilation operators 
denoted by b 1 ? and b 1 ^ are governed by the Heisen- 
berg equation with the diagonalized Hamiltonian H m = 

Eo(t) + ^ fe ^fe "fe- l n tnc diagonalized Hamilto- 
nian, the new ground state energy is defined by E (t) = 

R e (/fc W<Afc (*))• The energy spectrum hujk is 
that of the gapless excitations characterized by u>k — kc 
for a small fc, where c = gn/M is the speed of sound 
of the condensate. This is when k£ <C 1 holds, where 
£ = h/(2Mc) is the coherence length of the condensate. 
For large k, on the other hand, tkdk — £fe which is the free 
particle energy spectrum. We remark that the spectrum 
ujk and the speed of sound c are the same as in the origi- 
nal Bogoliubov model without impurities. Therefore the 
motions of impurity do not affect neither ujk nor c within 
our approximation. 

The vacuum state of the Fock space Hi n for b 1 ^ and b 1 ^ 
is defined by |0; n ) = 0. At remote past t = — oo, the 
vacuum |0i n ) coincides with the homogeneous condensate 
\BEC). The physical interpretation of the above result is 
quite simple. The motions of impurity creates the time- 
dependent coherent states. The dressed Bogoliubov's ex- 
citations do not annihilate the homogeneous condensate 
vacuum: 



b™ \BEC) = -4>k{t)\BEC). 



III. RADIATION SPECTRUM 



(15) 



One of main effects described by our model is the effect 
of energy radiation due to the motions of impurity in the 
homogeneous condensate. From our model we find that 
even though the systems is in the ground state where no 
elementary excitation is present, there is energy trans- 
ferred from the motions of impurity to massive bosons. 



In other words, the motions of impurity will dress bo- 
golons in such a way that bogolons will not see the ho- 
mogeneous condensate as the vacuum. Since we do not 
take into account any back reaction on the impurities, 
this transferred energy is the energy required in order for 
the impurities to move along given trajectories. It means 
that we need to feed this amount of energy to keep the 
impurities moving. More realistic model for the radiation 
effects in BEC should be done by considering quantum 
impurities and taking into account the dynamics of them. 

The occupation number for bogolons with respect to 
the homogeneous condensate 



h k (t) = (BEC\ b?b% \BEC) = 



n z \ z € k 
Nh 3 uj k 



h(t)\ 2 , (16) 



counting the emitted bogolons accompanying with the 
motions of impurity in BEC. 

We next evaluate the expectation value for the number 
operator Nk(t) — a k (t)ak(t) = a\.(t)ak(t) for the excited 
particles with respect to the vacuum (0i n |-/Vfc(i)|0j n ). This 
number counts the number of non-condensated particles. 
The depletion of the condensate d(t) due to the quantum 
fluctuation can be evaluated by summing over all modes 
and divided by the total particle number N: 



(17) 



+ 



N *-~> huk 

gn 



mt)-^- k (t)n as) 



Since \(/)k\ 2 has an additional factor 1/N, the motions of 
impurity do not contribute to the depletion of the ho- 
mogeneous condensate in the thermodynamic limit. In 
contrast, real experiments have always a finite number 
of particles and size. Therefore, there are effects on the 
depletion of condensates due to the motions of the impu- 
rity. 

The dissipated energy £k(t) for a given mode k is ob- 
tained by multiplying fife (t) by the excitation energy huik, 
i.e., £k(t) = hujkfik(t) and the total dissipated energy 
£(t) is given by summing over all modes. Therefore, for 
the t — > oo limit which takes into account the complete 
trajectory of impurity motion, we obtain rather a simple 
expression the total dissipated energy in the thermody- 
namic limit as 



lim V" £ k {t) 

t— >oo z — ' 

d 3 k 



h 2 



:C-k 



dt e iuJ kt-ikX(t) 



(19) 



When an impurity is moving with a constant velocity, 
i.e., — v ti the total dissipated energy per unit time 
is proportional to J d^ke^S^^ — k-v). This clearly shows 
that the speed of impurity v = \v\ needs to exceed the 
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speed of sound to create finite amount of radiation, oth- 
erwise there is no radiation from the impurity motion 
[2, 3, 6, 7]. When the impurity acquires some acceler- 
ation, on the other hand, the integral (19) produces a 
non-zero value, that is a finite amount of energy radia- 
tion. 



IV. RADIATION FROM ACCELERATED 
IMPURITIES 

We consider three examples to illustrate our model and 
its general solution. We are interested in the occupation 
number (16) and the total dissipated energy (19) for the 
t — > oo limit, i.e., in evaluating the integral: 



I k = lim I k (t) 

t^oo 



dt e 



iuikt— ife-C(t) 



(20) 



This limit, of course, should be understood as an idealiza- 
tion to obtain analytical results, and one needs to calcu- 
late a finite interval to compare with experiments. In the 
following, we also calculate the numbers of emitted bo- 
golons dh k between k and k + dk in the thermodynamic 
volume limit. Similarly, the energy radiated between k 
and k + dk is d£k = haj k dh k . We note that a constant 
translation to the impurity trajectory — »■ + Co 
gives rise to a phase factor which vanishes upon taking 
the absolute square. Thus any translation of the impurity 
trajectory will not produce any physical effect. 



A. Black-body-like radiation 

We first seek a trajectory which leads to the Planck- 
ian distribution for the occupation number (16). With 
a simple analysis, the trajectory with two parameters Co 
and T > 



C(t)=Coe- r °* (0,0,1) 



(21) 



leads to the integral 

/oo 
dt exp(iuj k t - ik z ( Q e- rot ). (22) 
-oo 

Here k z is the z component of the wave number vec- 
tor of the emitted bogolons. The parameter T can 
be expressed in terms of the acceleration constant a as 
To = a/c and we refer it to as the acceleration parameter. 
The case of negative acceleration T < can be analyzed 
similarly and we only consider the positive case in this 
paper. Changing the variable from t to r] = e~ r °* reads 



1 

r o Jo 



dvv -iw h /r -i e -ik M Cov m (23) 



To carry out this integral, one formally introduces the 
regularization uj k /TQ — > uj k /To + it where e is a small 



positive number which is set to zero at the end of cal- 
culations. Such a regularization, however, needs to be 
justified and this will be discussed in the next section. It 
is important to note that the z component of the wave 
vector k z can be both positive or negative. With the 
above mentioned regularization, we obtain the integral 



as 



J- o 



. i ^*) e -«n(fc.&)^ ) (24) 



where sgn(x) is the sign function and T(x) is the Gamma 
function. 

The appearance of the above integrals (22) is rather 
surprising when comparing with an informal derivation of 
the FDU effect using time-dependent Doppler shift, see in 
particular eq. (7) of Ref. [32]. Similar formulae have been 
used by many authors to discuss FDU-like effects in non- 
relativistic systems [9, 13, 14]. However our calculation 
exhibits a different result as discussed in the next section. 

The occupation number is 



n 2 X 2 e k 2tt 
NH 3 LJk WfcTi 



x 



*Jk Wfel 

l^/To.l]- 1 (fc,Co>0) 

{[l-e" 2 ™*/^]" 1 (fc,Co<0) 



(25) 



When the direction of emitted bogolon is on the xy 
plane, i.e., k z = 0, the above integral is proportional 
to the Dirac delta distribution 8(uj k ). Using the formula 
xS(x) — together with the factor e k /uj k , we get zero 
radiation for k z (o = 0. We now examine this result (25) 
for the infrared and the ultraviolet region in detail. From 
now on, we choose the impurity trajectory to be confined 
within the positive z space, i.e., Co > 0, and hence the 
impurity is accelerated from a positive infinity point +oo 
toward to the origin as the time changes from — oo to +oo. 
First of all, we observe that the distribution of emitted 
bogolons does not depend on the direction within the 
upper half volume (z > 0) and the lower half volume 
(z < 0), respectively. Secondly, the radiation is indepen- 
dent of a scaling with respect to the parameter of the 
trajectory Co and is determined solely from the accelera- 
tion parameter r . Lastly, the mathematical expression 
of the emission is the black-body radiation with the tem- 
perature ksTu = IiTo/2tt in the upper half volume and is 
completely different from thermal radiation in the lower 
half volume. We note a similar result was obtained and 
discussed in Ref. [13]. 

When the emission is restricted only for the low mo- 
menta, i.e., for the linear dispersion regime u> k ~ kc and 
kc <C r , the above two cases show the same asymptotic 
behavior. The number of emitted bogolons between k 
and k + dk in the thermodynamic volume limit is given 

by 



dn k 



X 2 



2ghc 



kdkdn for k < £ _1 ,r /c, 



(26) 
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where dfl is the element of solid angle. Although this for- 
mula resembles the Rayleigh-Einstcin- Jeans law for the 
black-body radiation, it cannot be interpreted as a sig- 
nature of the thermal radiation. The reason is simply 
because the effective "temperature" Tjj oc T docs not 
show up for the long wavelength approximation. 

A serious difficulty arises for the ultraviolet region 
where the emitted bogolons within an infinitesimal wave 
number dk is 



dfii 



^dunJ^f ^ >0) , 

h 3 T Q c 2 I i+ e -2»«*/ro (fc z £ < o) 



(27) 

for k 3> £ 1 , r /c. This becomes a constant and is in- 
dependent of the excitation energy of bogolons for the 
lower half volume as k approaches infinity. Therefore, 
the ultraviolet divergent for the total dissipated energy 
occurs for this half volume when integrating all over the 
wave numbers. 

One way to understand the origin of the above men- 
tioned ultraviolet divergence is to look at the form of 
potential corresponding to the trajectory (21). With a 
simple analysis of classical mechanics, the potential is 
V(C) = — MimpT^ 2 /2 with M imp the mass of impurity, 
which is an inverted harmonic potential. It is clear that 
the speed of impurity is not bounded and can become 
infinite for the asymptotic past and/or future. In reality, 
it is impossible to construct such a potential for the en- 
tire space since it requires an infinite amount of energy. 
Thus, this type of trajectories is only possible within a 
finite region and hence is limited by a given experimental 
situation, and one should not integrate the whole trajec- 
tory in the integral (20). 

For a finite integral interval specified with (£ i; tf), the 
integral of interest is calculated by the same regulariza- 
tion procedure as 



h(U, tf) 



\kz(o\ 



-sgn(fcz) TTU) k /2T 



7 H^fc z Co e - iot *)-7( 

1 



To 



ik z ( e- rot i) 



(28) 



where j(x, z) is the incomplete gamma function of the 
first kind. We note that in the limit — > — oo, the 
first term in (28) converges to the previous result (24). 
However, the second term oscillates with the angular fre- 
quency Uk as the final time tf approaches oo. This be- 
havior indicates that the regularization and infinite time 
integral limit do not commute each other and physics for 
finite time interval may be different from the above ide- 
alized infinite time interval limit. To analyze further, we 
consider a special case t% = and set tf — T. Using the 
asymptotic expansions of the incomplete gamma func- 
tion, the number of emitted bogolons takes in the long 
wavelength regime as 



dh k 



X 2 cT 2 
2gh 



k 3 dkdn for k < ^^q/c, (29) 



whereas it behaves as in the short wavelength regime 

dn k ~ 4n ^f c f [l-cosfc z Co(l - e" roT )] k~ 2 dkdfl, 

(30) 

for fc>r 1 ,r /c. 

Comparing with the previous expressions (26, 27), the 
finite interval case is totally different from the infinite 
case both for the long and the short wavelength regimes. 
Importantly, dhk does not decay exponentially for large 
value of wave numbers but decays algebraically. There- 
fore, the trajectory (21) cannot be used to observe the 
temperature-like radiation unless its acceleration covers 
an infinite time period. This is, however, impossible in 
real experiments and we conclude that the appearance 
of mathematical expression for the Planckian distribu- 
tion (25) is an artifact of invalid extrapolation for the 
physical parameters. 



B. Relativistic uniform acceleration 

The relativistic trajectory of a uniformly accelerated 
particle along the ^-direction is 



CW-^i + (f) 2 (0,0,1), 



(31) 



with a given proper acceleration a. To discuss a simu- 
lation of the relativistic effect in BECs, we replace the 
speed of light by the speed of sound c in the condensate. 
Note that this trajectory can be realized with a non- 
rclativistic potential of the form V(Q = ^Mj mp (c 3 /a() 2 . 
Importantly, the speed of the impurity is always bounded 
by that of sound, |C(i)| < c. In the limit c — > oo or within 
the short time period at/c <C 1, the trajectory becomes a 
non-relativistic quadratic one such as ai 2 /2+(constant). 
To evaluate the above integral (20), we change the inte- 
gration variable t to a dimensionless variable s through 
t = (c/a) sinhs, which is called a proper time in relativ- 
ity, and we define 

At fc = ~ a \pl~- (ck z ) 2 = ^\/ sin2 9 + fc2 £ 2 ' ( 32 ) 

(33) 

where 9 is an angle between k and the z axis. Then, 



<7fc = tanh X ( — ), 



2c 

h = — [7r £(/ife) cosh CTfc + zKi(/Zfc)sinh a*;], (34) 

a 

where Kj(x) is the jth modified Besscl function of the 
second kind. The first term in (34) will not contribute 
to the occupation number for the following reasons. The 
first is that ^ = if and only if k = 0. The second 
reason follows from two properties of the Dirac delta dis- 
tribution; S(f(x)) — '^2 i S(x — Xi)/\f'(xi)\, where Xi are 
the roots of f(x) — 0, and xS(x) = 0. The final reason 
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is the presence of the factor et/oJk in (16). Therefore the 
occupation number in the thermodynamic limit is 



2n\ 2 c 6 k 6 cos 2 6 (Ki{n k 
{2TT) 3 Ma A huj k 



dkdVl. 



(35) 



Let us now examine the result for the infrared and 
the ultraviolet region using the asymptotic behaviors of 
Kj (x) . The [ik becomes very small in the infrared region, 
and hence the leading term of (35) becomes 



dhk 



2nX 2 kcos 2 e 



(27r) 3 Mc 3 ft(sin 2 9 + fc 2 £ 2 ) 



-dkdfl. 



(36) 



Therefore, the leading term has a singularity around 
k ~ at the angles 9 — 0,ir, i.e., along the direction 
of trajectory, which results in the divergent result for 
the total number of emitted bogolons. This infrared sin- 
gularity due to bogolons has the same origin as the fa- 
mous infrared catastrophe in quantum electrodynamics 
[33]. There, the total number of emitted photons due to 
a motion of accelerated charge particle diverges as k —¥ 0. 
However, the physical measured quantity, the radiated 
energy, is still finite. 

We next use the asymptotic form of Ki for large /Xfe to 
get the leading term: 



dn k 



n\ 2 {2Mcfco& 2 9 
(2ir) 2 h 5 ak 2 



exp(- 



Ma 



)dkdn. (37) 



Hence the dissipated energy in the ultraviolet region after 
integrating over the range of solid angles is 



4n\ 2 M 2 c 3 . hk 2 c, „ 
d£ k ^ o_ t 3_ ex P( — Tl-) dk - 



Ma 



(38) 



Thus d£t t asymptotically depends on k through the ex- 
ponential factor exp[— hk 2 c/{Ma)\. We note that this 
factor can also be obtained by evaluating the integral us- 
ing the stationary phase method. This high momentum 
behavior could be used to measure the number hc/M. 

Lastly we notice that two length scales in this exam- 
ple have different orders of magnitudes; the coherence 
length £ and the acceleration length l a = c 2 /a. In most 
real experimental situations £ -C l a , since accelerations 
available in laboratories are very restricted. Then we 
again use the asymptotic form of Ki to estimate the to- 
tal dissipated energy for the weak acceleration limit as 



10ft 3 



n\ 2 M 2 c I ha 



irMc 3 ' 



(39) 



Although the speed of impurities never exceed the speed 
of sound we expect a small finite amount of radiation due 
to the uniformly accelerated impurity. 



V. DISCUSSION 

In this paper, we have investigated a simple model for 
the accelerated motions of classical impurity in the ho- 
mogeneous condensate and radiations accompanied by 



them. It is shown that the impurity couples linearly to 
the field of the Bogoliubov excitations as shown in (9). 
Since the effective theory of the low energy excitations 
are known to resemble a relativistic mass-less field the- 
ory [34] , one naively expects that this model can be used 
to simulate the FDU radiation for a suitably chosen accel- 
erated motion of impurity. Indeed, our model also shares 
similar mathematical formulas with other proposed mod- 
els as discussed in Refs. [9, 13-15, 32]. 

We have, however, pointed out that the appearance of 
the Planckian distribution for a radiation spectra corre- 
sponding to the impurity trajectory (21) is valid only in 
the mathematical limit in which the impurity is accel- 
erated from the infinitely past to the infinitely remote 
future. This limit requires infinite amount of energy and 
hence one can hardly expect it to be simulated in labo- 
ratories. To check whether or not any trace of such an 
extreme situation can be observed in real experiment, we 
have examined the finite acceleration time case in detail 
and have shown that the results are completely differ- 
ent from the idealized limit. Obtaining the Planckian 
distribution also requires an appropriate regularization, 
otherwise the integral does not converge to the expres- 
sion (24). This regularization can be understood, for 
example, as the adiabatic switching off of the impurity 
coupling, i.e., A — > Aexp(— et) with e — > at the end of 
calculation. We, however, emphasize that the usage of 
other adiabatic factors will result in different formulae 
and hence the appearance of the Planckian distribution 
is regularization dependent. Our result indicates that 
similar difficulties may exist for observing the Planckian 
distribution in other mentioned proposals. 

We have also examined a relativistic trajectory for a 
uniformly accelerated impurity and have analyzed the 
radiation spectrum in detail. The result (38) shows an 
exponential decay of the emitted excitations for short 
wavelength behavior of the radiation spectrum. This ex- 
ponential factor is sometimes referred to as a Boltzmann 
factor of the Bose-Einstein distribution and is related 
to the effective temperature of the FDU effect. With 
this correspondence, the effective temperature here is 
ksT — ha/2c which differs by a factor it from the orig- 
inal FDU temperature Tjj = hT^/2iT. It is questionable 
if we can relate an appearance of the Boltzmann factor 
to the black-body radiation which is an incoherent radia- 
tion from the thermal equilibrium state. This is because 
radiations from the motions of impurity in condensate 
are coherent processes. The same statement also holds 
for the FDU effect. The thermal nature of reservoir is a 
consequence of a statistical ensemble average. From this 
averaging procedure, the expectation value of occupation 
number for bosons results in the Planckian distribution, 
but not the other way around. In other words, an ap- 
parent mathematical expression does not mean that the 
system is in a thermal equilibrium at certain tempera- 
ture. This point was originally discussed in Ref. [35]. 

Along the same line of analysis, we gave a preliminary 
study for a radiation spectrum from a constant circular 
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motion [7, 8]. It is expected that this kind of impurity 
motion can be realized in the current experimental tech- 
niques. In future, we intend to study the circular motion 
in more detail based on realistic settings such as including 
a trapping potential, finite temperature effects, the back 
reaction to the impurity, and so on. It is also interest- 
ing to examine a connection to the so-called the circular 
FDU radiation and to discuss a possibility to simulate 
such an effect in BECs. 

Finally, we briefly comment on a possible connection 
to the emergent theories of gravity. Several authors have 
proposed that gravitation is considered as an emergent 
phenomenon arising from collective excitations in quan- 
tum many-body systems [36-45]. In these proposals a 
physical vacuum is considered as a real condensed state 
which contains a huge number of heavy bosons, and the 
Minkowski vacuum is assumed to be the homogeneous 
condensate. It is clear that effective theory breaks down 
at a certain energy scale and the deviations from the 
prediction of the relativistic quantum field theory are ex- 
pected to occur. In this line of research, one is interested 
in these deviations and the new phenomena predicted 
by the microscopic theories behind, that are possibly ob- 
served in astrophysics. 



We are then interested in looking at the radiation for 
beyond the linear dispersion regime. In the relativistic 
uniform acceleration, for example, we see that the short 
wavelength behavior of the radiation spectrum is propor- 
tional to exp(— fik) — exp(— cjfec/a) along the direction of 
the trajectory. This dependence can be used to estimate 
higher modes of the excitation spectrum of the Bogoli- 
ubov excitations. It is also an open question to investi- 
gate if one can simulate the FDU effect itself in BECs, in 
which any accelerated observer perceives the Minkowski 
vacuum as a thermal state, and to clarify the physical 
mechanism leading to this effect if it exists. 
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